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Abstract. In this paper the author studies the problem of the homogeniza- 
tion of a diffusion perturbed by a periodic reflection invariant vector field. The 
vector field is assumed to have fixed direction but varying amplitude. The ex- 
istence of a homogenized limit is proven and formulas for the effective diffusion 
constant are given. In dimension d = 1 the effective diffusion constant is al- 
ways less than the constant for the pure diffusion. In d > 1 this property no 
longer holds in general. 



We consider the problem of the homogenization of a diffusion perturbed by a 
reflection invariant vector field. The general set up we have in mind is to understand 
the limit as e — > of the solutions u e to an elliptic equation, 

(1.1) — — A u e {x,oS) — 2b £ (x,uj)du e (x,oj)/dxi 



In l|l.l(l the function / : R d — ► R is smooth of compact support and SI is a probabil- 
ity space. For simplicity we have assumed that the vector field is always in the x± 
direction and hence can be described by the scalar function b e . As e — > the field 
becomes rapidly oscillatory and therefore one might expect that u £ {x,oj) converges 
with probability 1 as e — > to a homogenized limit u(x) which is the solution to a 
constant coefficient elliptic equation, 



The effect of the rapidly oscillating vector field b £ is contained in the coefficient q(b) 
in lfTT2)l . 

In order for a limit u{x) satisfying H1.2(l to exist it is necessary to make as- 
sumptions concerning the rapidly oscillating field b e . These are primarily that 
the distribution functions of the variables b e (x, ■), x £ R d , are translation and re- 
flection invariant. To be specific, we assume that there are translation operators 
t x : $1 — > SI, x £ R d , which are measure preserving and satisfy the group prop- 
erties T x T y = r x+yi x,y £ R d , t ^identity. Suppose b : SI — > R is a bounded 
function. We then set b e {x,u) = b{r x / e to), x £ R d , u> £ fi, e > 0. Such a b e has 
translation invariant distribution functions and is rapidly oscillating as e — * 0. For 
b £ to satisfy reflection invariance we let R : R d — > R d be the reflection operator 



1. Introduction 



2d 



+ u E (x,u) = f{x), x = (xi, ...,x d ) £ R , u> £ SI. 



(1.2) 
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R(xi, Xd) = (—xi,X2, Xd), x = (xi, Xd) € R rf - Wc then require b : — > R 
to satisfy the identities, 

(1.3) ^n&fo* •)) = (n^* y > €Rd ' 1 - n > ™- 

where (•) denotes expectation on f2. Evidently (|1.3|) implies that (&(•)) = 0, so the 
vector field has no net drift. 

A concrete example of an fi and a 6 : £1 — * R which satisfies JO} is given by 
taking to be a torus, 17 = Y\ i=1 [0, Lj] with periodic boundary conditions and 
uniform measure. The operators t x : Vl — > f2, x € R d , are just translation on il and 
reflection invariance of (|1.3[1 is guaranteed by the condition, 

(1.4) 6(xi,x 2 , ...,x d ) = -b(Li - xi,x 2 , ...,x d ), x = (xi, x d ) S fi. 

We shall show that for a discrete version of a periodic fl with 6 satisfying l|1.4|l a 
homogenized limit exists with q(b) satisfying < q(b) < oo. For d — 1 one has 
#) < 9(0) = 1/2. For d > 1 it is no longer the case that q(b) < q(0) = l/2d 
in general although this does hold for L\ sufficiently small. One might wish to 
understand this difference between d = 1 and d > 1 by observing that only in d > 1 
can one construct nontrivial divergence free vector fields. The homogenized limit of 
diffusion perturbed by a divergence free vector field necessarily yields an effective 
diffusion constant which is larger than the constant for the pure diffusion [Sj. 

The homogenization problem considered here appears to have only been studied 
in the case where Q is an infinite space for which the variables b(r x ■), x S R, are 
uncorrelated on a scale larger than O(l). The problem was introduced by Sinai 
[T7j in a discrete setting. He proved that in dimension d = 1 a scaling limit of the 
random walk corresponding to a finite difference approximation to Ijl.f I) exists with 
probability 1 in SI. The limiting process is strongly subdiffusive. In a subsequent 
paper Kesten obtained an explicit formula for the distribution of the scaling 
limit. For dimension d > 3 Fisher ^U] and Derrida-Liick [5] predicted that a 
homogenized limit exists as in l|1.2fl with < q(b) < oo. This was proved for 
sufficiently small b by Bricmont-Kupiainen |S] and Sznitman-Zeitouni |2U) using a 
very difficult induction argument. A formal perturbation expansion for q(b) was 
obtained in [S] where it was shown that each term of the expansion is finite if 
d > 3. One does not expect the series to converge however. For d — 1,2 there are 
individual terms in the perturbation expansion which diverge. 

A main difficulty in the homogenization problem (|l.l(l . (jl.2|l is that when Q, is 
infinite, good a-priori estimates on the solution to ll.lfl do not hold for all config- 
urations of &(•). In contrast such estimates do hold for divergence form equations 
with zero drift. The proof of homogenization in these cases is therefore consider- 
ably simpler than for the problem (|1.1|) , i|1.2fl • The first proofs of homogenization 
for divergence form equations were obtained by Kozlov 12; and Papanicolaou- 
Varadhan £ 15; in the continuous case. Kiinneman |13j proved a corresponding 
result for the discrete case. For non-divergence form equations with zero drift the 
first proofs in the continuous case were given by Papanicolaou- Varadhan |16| and 
Zhikov-Sirazhudinov [221 ■ Lawler and Anshelevich et al pQ proved homogeniza- 
tion for a discrete version. See the books of Bolthausen-Sznitman [2] for an account 
of the theory in a discrete setting and of Zhikov et al [23 f° r the continuous case. 
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In this paper we shall be concerned with a discrete version of the homogenization 
problem described by Ijl.ljl . Q1-2JI . (|1 .3|> . Thus the probability space 17 is acted upon 
by translation operators t x : 17 — > 17 where now x £ Z d , the integer lattice in R d , 
and satisfy the group properties t x r y = r x+y , tq = identity. For i — l,...,d let 
ej £ Z d be the element with entry 1 in the ith position and in the other positions, 
the discrete equation corresponding to is given by 

d 1 

(1.5) u e {x,Lo) - 2_j ~2d i u e( x + zeiiu) + u £ (x - eei,uj)] 

i=l 

- K T x/e) [u-s(x + sei,w) + u £ (x - eei, u)] 

+ e 2 u e (x, lu) = e 2 f{x), x e Z d = e Z d , uefl. 

We assume that b : 17 — > R satisfies sup w \b(uj)\ < l/2d, in which case ()1.5|l is an 
equation for the expectation value of a function of an asymmetric random walk. 
Hence ()1.5f) has a unique bounded solution. We assume that b satisfies the reflection 
invariant condition (|1.3(l (with Xi S Z d , 1 < i < n, now). We also assume that il is 
finite, in which case one can see (Lemma 2.4) that 17 is isomorphic to the integer 
points on a torus and b has the reflection invariance property l|1.4|) . In §2 we prove 
the following theorem (with |_-J denoting the integer part): 

Theorem 1.1. Assume fl is a finite probability space and the translation operators 
t x : n — > n are ergodic, x £ Z d . Then there exists q(b), < q(b) < oo such that 
with u s the solution to \1.5)) and u the solution to H.ty) . 

lim sup \u £ (e[x/e\,u}) ~ u(x)\ = 0. 

Suppose now that 17 consists of the integer points on the torus ]X=i ^i] C R d 
with periodic boundary conditions. The reflection invariant condition correspond 
to (EH is given by 

(1.6) b(xi,x 2 , .:,x d ) = -b(Li - 1 - xi,x 2 , —,x d ), x = (xi,x 2 , ■■■,x d ) £ 17. 

For b satisfying (|l.t)|) we prove in §2, §3 the following results concerning the coeffi- 
cient q(b) of the homogenized equation (|1.2|l : 

Theorem 1.2. (a) For d = 1 one has q(b) < 1/2. 

(b) Ifd>\ and Li =2 one has q(b) < l/2d. 

(c) If d — 2 and L\ — 4 one has q(b) < 1/4. 

(d) If d > 2 and L\ > 6 is even then there exists b with q(b) > l/2d. 

The proofs of (a), (b), (c) are given in §3 and are based on applications of the 
Schwarz inequality. The proof of (c) is quite lengthy and depends crucially on actual 
numerical values for a Green's function associated with standard random walk on 
the integers. The proof of (d) is given in §2. One observes that perturbation theory 
yields q(b) = l/2d + 0(|6| 2 ) and that the term 0(|6| 2 ) can be positive. 

In the proof of Theorem 1.2 we use a representation for q(b) in terms of invariant 
measures for random walk on 17 with drift b. Let £ld-i consist of the integer 
points on the d — 1 dimensional torus 11^=2 with periodic boundary 

conditions. Setting L\ = 2L with L an integer we define 12 by 

Q = {(n, y) : 1 < n < L, y£ n d -i} , 
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whence fi is the double of f2. Observe that the boundary of dtt is given by 

dh = {(l,y), (L,y):y efid-i}. 

Let ip* be the invariant measure for random walk on Cl with drift &(•) in the ei 
direction and reflecting boundary conditions on dCl. We define ijj : Qd-i — > by 

= [l/2d- 6(1, y)] ¥>*(l,y), y e Q d -i. 

Then q(b) is given by the formula, 

(1-7) q(b) = U{-A d _ 1 +4]- 1 ^ R ) 

where V'i? is defined exactly as tp but with 6 replaced by —6. In l|1.7|l the expectation 
(■)q i is the uniform measure on fid-i and A^-i is the d — 1 dimensional finite 
difference Laplacian on functions with domain f^-i- The normalization of tp is 
chosen so that q(0) = l/2d. The general formula ljl.7|l is proven in §4. 

2. Proof of Theorem 1.1 

We follow the method introduced in to obtain homogenized limits. Thus in 
(11.5(1 we put u £ (x,uj) — v E (x,T x / e uj) whence 1)1 .5|l becomes 

2d ^ 

(2.1) v s (x,uj) - ^ — [v e {x + eej,T ez uj) + v £ {x - £e t ,T_ e , uj)] 

i—l 

- b(uj) [v £ (x + eei , r ei uj) - v e {x - eei, r_ ei w)] 

+ £ 2 « e (a?,w) = £ 2 /(«), xeZ d = eZ d , uj e (1. 

Next we wish to take the Fourier transform of (|2.1|l . To show that this is legitimate 
we first show that the solution u £ (x,uj) of l|1.5fl decreases exponentially as x — > oo. 

Lemma 2.1. Suppose f : Z d — > R /ias finite support in the set {x — {x\ 1 ...,Xd) £ 
Z d : \x\ < R}. Let u £ (x,uj) be a bounded solution to Then there are constants 

C, K(e) > such that 

(2.2) \u £ (x,uj)\<CeMK(£)(R-\x\)]\\f\\oo, xGZ d . 



Proof. We write u £ (x,uj) 
satisfies 



~ kxi u £ ,k(x,uj). Then from l|1.5|l the function 



U £ ,k 



(2.3) 



1 d 

— ^2 I 2u £:fc (x, uj) - u £tk (x + eei , uj) - u £tk (x - eei , uj)] / e 2 



-ke 



Td +b(r x/£ uj) 



{u £:k {x,uj) ~ u £tk {x + eei,uj)} /e 2 



ke 



Yd~ h{T ^ u) 



{u £ik (x, uj) - u £ . k {x - eei, uj)} / e 2 

+ {l — [coshfe — l]/efe 2 + 2b{r x / £ uj) sinh ke/e 2 } u £t k(x, uj) 

= e kxi f(x), x € Z d . 

We may assume wlog that / is nonnegative, whence u £ ^k is also nonnegative. Sup- 
pose u £ ^k attains its maximum at a point x £ Z d . Then we have that 

(2.4) {l - [cosh ke - l]/de 2 + 2h(r x/e uj) sinh ke/e 2 } u Sik (x, uj) 

<exp[k(x • ei)] H/lloo, \x]<R, 
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whence it follows that 

(2.5) u e (x,tu)<Cexp[k(x-x)-e 1 ]\\f\\ 00 , x G Z d . 

We need to show that the point x exists for sufficiently small k. To see this assume 
for contradiction that it does not exist. Then l|2.3[) implies that snp^ x ^ <N u £ ^(x,uj) 
grows exponentially in N as N — » oo. The rate of exponential growth remains 
bounded away from as k — > 0. Hence, taking k sufficiently small, we conclude that 
the function u e is unbounded, contradicting our assumption on u e . The inequality 
(I2.2|l now follows from l|2.4(l , (|2.5|) on generalizing to all directions ey , 1 < j < d. □ 

For £ G [—tt/e, Tr/e] d we put 

(2.6) « B (C,w) = f v e (x,oj)e ix <dx= ^ e d v £ (x,cu)e lx < . 

Then from l|2.1|l we have that 
d 

(2.7) 6 e (£,a;) - ^-[e-^z) e (e,Te^) + e ^« £ (e,T_ e ^ 

j=i 

- b{cj) [e- is ^v e (Z,T ei u) - e K ^ E (C,r_ eiW )] + e 2 v e (Z,Lj) 

= e 2 f £ (0, £ G [-tt/e, 7r/e] d , w6 S], 

where f e denotes the discrete Fourier transform l|2.6|) of /. To solve l|2.7l) we define 
for C G [— 7r, 7r] d an operator £^ on functions ^ : £1 — > C defined by 

d 

(2.8) £ c *(w) - [ e -^*(r e .a;) + e^*(r_e J .a;)] 

3=1 

Next we define an operator T v ^, 77 > 0, ( G [— vr,7r] <i on L°°(f2) by 

(2.9) T Vt £ ip(uj) = 77 [£f + rj\ 1 ip(oj), we!!. 

It is easy to see that T Vi c, is a bounded operator on L°°(f2) with norm at most 1. 
In fact the RHS of l|2.9[) is the expectation for a continuous time random walk on 
il x Z d . The walk is defined as follows: 

(a) The waiting time at (00, x) G f2 x Z d is exponential with parameter 1. 

(b) For j = 2, d the particle jumps from (lj, x) to {T ej oj, x + ej) with proba- 
bility l/2d and to (j— e ,u),x — e^) with probability l/2d. 

(c) The particle jumps from (u>, x) to (r ei u>, x+ei) with probability l/2d+b(u>), 
and to (r_ ei oj, x — ei) with probability l/2ci — 6(w). 

If [u)(t),X(i)] G f2 x Z d is the position of the walk at time t then 
(2.10) 

'POO 

T rh( ip(uj) =rjE dt e^VM*)) exp[-iX(t) • C] w(0) = w, X(0) = . 

It is clear from the representation i|2.10|l that ||T^ fHoo < 1. We conclude from this 
that l|2.7|) is solvable with solution given by 

(2.11) e e ({,w) = f s (0 T^ M (l)(u), u G n, g G [-7r/e,7r/ £ ] d . 
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To obtain the homogenization theorem wc need then to obtain the limit of the RHS 
of l|2.11(l as e — > 0. To facilitate this we observe from Ij2.8(l that 



1 

(2.12) [£ c +?7]1 = 1 + 77- -^cosCj + 2i6(w)sinCi. 

It follows therefore that 



(2.13) 



1 

1 + 7?- COS Cj 



3=1 



2i sin Ci 



1 d 

1 + T '- 3E cos ^ 



j'=i 



T, ; , c 6(w). 



Setting i] — e 2 , ( = e £ for some fixed £ € R d we see from (I2.13|) that 



(2.14) 



hmT £2 ^(l)H = I 



2^i 



1 rf 

1 + 5E« 



lime 1 T e 2 e j6(a;). 



We shall show that under the assumption of (jl.3JI the limit on the RHS of (|2.14(1 
exists. To do this we define two subspaces of the space L°°(n). We define L^(fi) 
as all functions $ S L°°(Q) such that 

U(r x -)p^b(T Xi -)\ = (-1)" +1 ^HT Rx -)f[b(T Rx ,)^ , aj,Si € Z d , 1 < i < n,n = 0, 1,2, ... 

Evidently (JOJ implies that 6 S L^(Q). We also see that if $ e Lg> (ft) then $(r ej •) 
and <&(r_ ej -) are also in L^(fi), j = 2, For j = 1 one has that $ £ L^(Q) 

implies both [$(r ei -) + $(r_ ei 0] and 6(-)[$(r ei -) - $(r_ ei ')] are in if(O). The 
space L^(f2) is defined similarly as all functions $ S L°°(f2) such that 

/^xOn^O) = ('IT (^(TRx^pbiT^J ,X, Xi E Z d ,l <Z < 71, 71 = 0,1, 2,.... 

From Q1.3[l we see that the function <E> = 1 is in L^(fl). As for the space L^(£l), 
if $ £ L^(fl) then $(r ej -) and $(T_ ej -) are also in £^(0), j = 2, Similarly 
both [$(r ei -) + $(r_ ei -)] and &(-)[$(r ei -) - $(r_ ei -)] are in L%(to). We note that 
the mapping <&(•) — > 6(-)<I>(-) maps i^(il) into 1/^(0) and vice-versa. 

We denote the operator of l|2.8|l for ( = by C. It is evident that C is 
the generator of a random walk on fl. Hence the kernel of the operator C is just 
the constant function. Furthermore C leaves the space L^(ri) invariant. Since 
the constant function is not in L^(O) it follows that there is a unique function 
ip € L^(n) such that 

(2.15) Cip = b. 

Let tp* be the invariant measure for the walk on f2 generated by C. Thus ip* > 0, 

(2.16) £V = 0, (<P*) = l, 
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where C* is the adjoint of C. Since C is non singular on the space L^(S1) it follows 
that tp* is orthogonal to L^(O). We can also see that tp* G L^(il). One simply 
notes that both C and C* leave the space L^(fl) invariant and that the constant 
function is in L^(S1). We obtain the limit on the RHS of (|2.14|l in terms of the 
functions tp, tp* defined by (|2.15l) , 12.1611 . 

Lemma 2.2. Let E L°°(tt) be defined by 

(2-17) V(-) = { 7^ + KO 1 <PM b (-)\ ^-r), 



^2d K ' } ry 1 ' [2d 
where tp is given by \2.15)) . Then if tp* is as in \2.1b)) there is the limit, 

d 

(2.18) lim e^T^^u) = < tp*i> > j 



2d^ Sj 

for all u> G f2, provided £i is sufficiently small. 

Proof. For ry > 0, £ £ [— 7T, 7r] d , let (^(ry, £) be the unique solution to the equation 

(2.19) [Cc + v]<P(v,Q =6- 
It is clear then that 

(2.20) e- 1 T E 2^b = £tp(e 2 ,eO- 

We define operators Af, by = + Cr^]/2, — [Cq — Crq]/2, where 
i?(Cij (d) — (— Ci> C2> •••) CO- The operator leaves the spaces L]f (ft) , L"^ (ft) 
invariant whereas B^ takes L^(fl) into L^(f2) and vice versa. Equation l|2.19|l is 
now equivalent to 

(2.21) [(A c +r ] ) + B c ]p(r 1 ,0=b , 

and we may write the solution of this formally as a power series, 

00 

(2.22) tpin, C)=E {~( A < + Vr'BcV (A c + V)- 1 b. 

The operators B^, (A^ + n)^ 1 on L°°(0) have norms satisfying \\Bq\\ < C2IC1L 
\\(Ac + if) \\ < I/ 7 ?; f° r some constant C 2 - Since Ao = £ is invertible on L^(U) it 
follows that for (77, ^) sufficiently small the operator norm of (Aq + n)^ 1 acting on 
L^{yi) satisfies + < Ci for some constant Ci. We conclude therefore 

for (n, () sufficiently small that 

(2.23) {(Ac + ny'BcYiAc+ny'b < C^|Ci| n C^^pU , 

00 

where r = n/2 if n is even, r — (n + l)/2 if n is odd. Hence if (77, C) an d ICil 2 / 7 ? 
are small then the series in (I2.22|) converges in L°°(£l) to the solution of l|2.21|l . 
It follows that for £ £ R d fixed with £1 sufficiently small we may construct the 
function tp(e 2 ,e£,) by means of l|2.22|l as e — > 0. 
To obtain the limit in (|2.18l) we write 

(2-24) tp(r),() =tpx{v,0+Mv,0 

where tp\ (77, Q) is the sum on the RHS of l|2.22|l over odd powers of n. It is evident 
from $F1% that for |£i| < l/C 2 \ft\ one has 

(2.25) lim etp 2 (e 2 ,e£) = 0. 
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We consider the first term in the sum for ip-y. Setting 77 = e 2 , Q = e£ and multiplying 
the term by e as in (|2.2Q(I we see that 



lim e (A ei + e 2 ) 1 B ei + e 2 ) 1 b 



lim i£ 2 £,\ (A e £ + e 



2\-i 



2rf 



2d 



where cp is the solution of 1)2. 15[) . Observe now that for any rp £ L°°(f2) we have 

- limie 2 £i (A e t + e 2 )~ 1 -(/> = -i£i(<p*i/j) lime 2 (A e < + e 2 )~ 1 1 , 

where 99* is the solution of l|2.16[l . From 12.12|l we have that 

d 



(M+e 2 )- 1 l=e 2 



l+e z 



3=1 



whence we conclude that 

lime 2 (A e£ +e 2 ) _1 1 = 1 



— Ve 2 



t2 

2d^-" J 

3=1 



We have therefore obtained a formula for the limit as e — > of the first term in 
the series representation of eipi (e 2 , e£) . Using the same argument we can obtain a 
formula for all the terms. For the rth term corresponding to r = (n + l)/2 with n 
as in (|2.22|) we see that the limit is given by the formula, 



(2.26) 



1 d 

3=1 



where -0 is the function (|2.17|l . Evidently if> £ L ( S(Cl). We have already observed 
that (f* is also in L^(S7). We conclude that 



lime<pi(e 2 ,eO = -i^(cp*ip) j 



3 = 1 



Then (|2~T%|) follows from this and (|2~2^|) . □ 

Lemma 2.1 enables us to compute the limit (|2.14f) when £1 is small. We have 



(2.27) 



hm7>^(l)H = 1 



6^0 



2d 



3=1 



We wish now to extend the identity l|2.27jl to all £ g R d . 

Lemma 2.3. Let K C R d be a compact set. Then the limit is uniform for 

£ e K, cj e Q. 

Proof. Since the LHS of I2.27fl docs not exceed 1 in absolute value we conclude 
that 



(2.28) 



- + 2(^>>0. 
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The inequality (|2.28(l in turn implies that the expression Ij2.26|l is the rth power of 
a number strictly less than 1 provided we also assume that 

(2.29) 2 (<p*tp) < l/2d. 

We show that the power series methods of Lemma 2.1 apply to prove the result 
under the additional assumption l|2.29[l . We shall see in §3 that Up*^) < for 
dimension d = 1, in which case (|2.29|) certainly holds. Since the constant function 
is the unique eigenvector of Aq = C with eigenvalue and it is also an eigenvector 
of Aq it follows that there exists S > such that if |C| < 8 then the adjoint Aq of 
Aq has a unique eigenvector tp*^ with eigenvalue equal to the eigenvalue of Aq for 
the constant function. Normalizing ip^ so that < ip^ >= 1, it is easy to see that 
there is a constant C\ such that 

(2.30) IK-p1°o<CiICI, ICI<* 

For |£| < 5, r\ > we define a projection P nt Q by 

p vX </> = (<p\ 1>) (Aq + t?)- 1 i, v e L°°(Cl). 

Then there is a constant C 2 such that 

(2-31) wiAQ + nr'-p^w <C7 2 , |C|<& 

The uniform convergence of (|2.18() for £ E K follows now 

from JOB, lESO just as 

in Lemma 2.1. 

Finally we consider the situation where l|2.29|) is violated. As in Lemma 2.1 we 
decompose the solution <p(r]X) of H2.19(l into a sum H2.24j) . The function <p>\{q, C) 
is a solution to the equation, 

(2-32) [Aq + V -Bq (Aq + ry)- 1 Bq] Vl ( V , Q = -Bq(Aq + i^ 1 b. 

The function ^2(77, C) is a solution to the equation, 

(2.33) [Aq + V -Bq (Aq + Ty)- 1 Bq\ <p 2 ( v , C) = b. 

It is easy to see that if <p2(r], £) is a solution of H2.33fl then the function (pi(rj, () = 
viViC) — ^2(77, C)i where <p(r), £) solves (I2.19|) . is a solution to l|2.32|l . Hence if 
(|2.32ll . (|2.33|l have unique solutions <p>i(r], f), (p2(r), C) then the identity l|2.24|l holds. 

We show that fTSfy has a unique solution in (J7) provided rj > and (77, £) 
are sufficiently small. To see this we write (|2.33|l as 

(2-34) [Aq + V - L rhC - B ( P vX Bq] <p 2 ( V , C) = b, 

where by (|2.31|) the operator L vX is invariant on L^(f2) and satisfies ||£j7,f || < C | C 1 2 
for some constant C. Next let (^3(77, C) be the solution to 

(2.35) [A ( + r l -L ri ,Q}<p 3 (r],0 = b. 

For (77, C) small there is a unique solution to (|2.35(1 in L^(fl) which satisfies 
(2-36) y~Mv,0\\oo <C[\r 1 \ + ICI 2 ], 

for some constant C, where <p is the solution of (|2.15|) . Now it is easy to see that 
the solution tp 2 (r], C) of l|2.34|) is given in terms of c/?3 (77, £) by the formula, 



(2.37) ^,(77,0 = 



1 * 



^3(7?,C), 
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1 d 

1 + r l- -j^2 C0S <j + 2 (<Pc B <Pa(V,Q) sin 2 Ci 



i=i 

where the operator B is defined by Bq = ism(iB. In view of (|2.36|) and the fact 
that Bip = ip and we are assuming (|2.29[) is violated, it follows that the denominator 
in (|2.37|) is positive for (77, Q sufficiently small. We have shown a solution ^2(11, C) of 
(|2.34|l exists in (f2) . The uniqueness of the solution follows from the uniqueness 
of the solution to (j2~35|) . Evidently the limit lj2~2"5|) follows from (JOB, fFTfy for 
all £ and is uniform for £ restricted to a compact subset of R d . 

Next we show that <|2.32[) has a unique solution in L^(f2) provided 77 > and 
(rj,£) are sufficiently small. First note that for (77, £) small the operator Bq(A^ + 
7])~ 1 Bq leaves L^?(f2) invariant and there is a constant C such that 



(2.38) 



We define the subspace £q of Lg'(fl) by 

Let be the projection operator on L^(r2) orthogonal to whence 
Consider now the equation related to (|2.32[1 given by 



(2.39) 



<P4.(V>0 



A c +7 1 - (I-P c )B c (A c +r))-'B ( 

= -(i-p^iAc + nyH. 

In view of l|2.38|l it is clear that for (77, sufficiently small the equation 
a unique solution (^4(77, £) in Furthermore, if we define </?i(t7, £) by 



has 



(2.40) Mv,0 



1 

i=i 



^ 4 (r/,C)-*sinCi(^S(A c +77) 
- sin 2 Ci (<p* c B(Ac + v)- 1 B ^4(77, 0) } / 

r 1 d 

1 + ^ - 2 ^ cos Cj + 2 sin2 Cl ^c B (A + 7/T 1 6) 

- 2z sin 3 Ci (<P* ( B ( A C + 77)- 1 B<p 4 ) 



then one sees that the formula (|2.4U|I yields a solution to H2.32[l . Conversely, since 
we are assuming (|2.29|l is violated, it follows that for (77, £) small (|2.4UI) is the 
unique solution in L^(n) to (|2.32|) . It is easy to see now from (|2.4U|) that the 
limit lim e ^o eipi(e 2 , e^) exists and is uniform for £ in a compact subset of R d . 
Furthermore, the limit is given by the RHS of (|2.18(1 . 

Finally we show that if <pi(rj, Q,ipv{ri, Q are solutions to lj2.32|) . I|2.33|l then (|2.24(l 
holds. To see this we put <p(r],£) = fiivX) +<P2 (?7, C) an d note that (|2.32|l . (|2.33|l 
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imply that ip(j], C) satisfies the equation 

[A c + r, - B C (A C + vr'Bc] tpfa Q=b- B C (A C + rj)- x b. 
We can rewrite this equation as 

[A c + B C + V - B ( (A C + ry)" 1 ^ + B c + v )] tp( v , ()=b- B C (A C + n)-\ 
which is the same as 

[C RC + V ] (A ( + r ] )- 1 {£ ( +T 1 }<p(i 1 ,0 = [C RC + V] i A C + v)' 1 b. 
Now using the fact that the operator Cr^ + rj is invertible we obtain (|2.19|) . 



□ 



Next we show that there is strict inequality in l|2.28[l . In order to do this we 
shall first obtain a concrete representation of spaces f2 which satisfy i|1.2fl . 

Lemma 2.4. Let Q, be a finite probability space and b : fl — > R satisfy \1. 'jj) . 
Then f2 may be identified with a rectangle in Z d with periodic boundary conditions. 
The operators t x , x £ Z d , act on f2 by translation and the measure (■) is simple 
averaging. Let R : CI — > £1 be the reflection operator defined as reflection in the 
hyperplane through the center of Q with normal e±. Then there is the identity 
b(uj) = -b(Rw), wefi. 

Proof. Since SI has no nontrivial invariant subsets under the action of the T ej , 1 < 
j < d, it is isomorphic to a rectangle in Z d with periodic boundary conditions. 
Thus we may assume fl is given by 

(2.41) n={x=(xi, 

where L\, ...jLd are positive integers. The action of the r ej is translation, r ej x 
x + Bj with periodic boundary conditions. The measure on fl is averaging, 

(2.42) (*(.)) = _J__ Yl •»,*«*)• 

' '' ■ :xi<Li-i, 

Ki<d 



, x d ) £ Z d : < Xi < Li - 1, 1 < i < d} 



L\Li 



Functions \& : £1 — > R are isomorphic to periodic functions if? : Z d — ► R. 

Next we consider the condition 11.3fl . We define a function b R ■ : £1 — ► R 
by b R (uj) = —b(Ruj), uj £ 0. It is easy to see that b R {r x uj) = —b{Rr x uj) — 
—b(TR X Rui),ui £ £1. Since R leaves the measure (|2.42f) invariant (|1.3|) implies that 
for any 9\,...,9k £ R, x%, ...,Xk £ Z d , there is the identity 



exp 



k 

E^ 

i=i 



Kr X] -) 



= ( exp 



k 

E 



We conclude that b= b R . 



□ 



Next we wish to construct the solutions tp, ip* of H2.15|l . (|2.16|) on the domain 
£1 defined by l|2.41|l . First observe that since £1 is the fundamental region for 
the homogenization problem we can assume that L\ is an even integer by simply 
doubling £1 if Li is odd. In that case the function b is determined by its values 
b(x), x £ CI, < x\ < L\/2 — 1. Hence we define a new fundamental region £1 by 



(2.43) 



fi = {x £ fi : < xi < Li/2 - 1}. 
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Wc can extend functions $ : f2 — > R to f2 by either symmetric or antisymmetric 
extension. For a symmetric extension we define ^ on fl — Q by 

(2.44) V(xi,...,x d ) = V(Li - 1 - xi,x 2 ,...,x d ), L 1 /2<x 1 <L 1 -l. 
For an antisymmetric extension we define \t by 

(2.45) $>(xi,...,x d ) = - f -xi,x 2 ,...,Xd), Li/2 < xi < L x - 1. 

Lemma 2.5. T/ie solution ip : £1 — > R o/ is an antisymmetric extension of 

its restriction to ft. The solution ip* : £1 — > R o/ 1^2. lb}) is a symmetric extension 
of its restriction to f2. 

Proof. This follows easily from the fact that 6 : f2 — > R is an antisymmetric exten- 
sion of its restriction to tl and the uniqueness of the solution to l|2.f 5J1 , l|2.16(l . □ 

Lemma 2.4 implies that we can find the functions (p, ip* by solving 12. 151) . (|2.1ti|) 
on fl with antisymmetric and symmetric boundary conditions respectively. Thus C 
acting on functions \f" : O — * R with antisymmetric boundary conditions is defined 

by 

d 1 

(2.46) £*(x) = *(x) - X/ 2d + ^ + *^ ~ ^ 

i=i 

- 6(x) [$(a; + ei) - &(x - e^] , x E £l, 
where the boundary conditions are given by, 
(2.47) 

^(-l,x 2 ,...,x d ) = -y(0,X2,...,x d ), V(Li/2,X2,...,x d ) = -$(Li/2-l,X2,...,x d ), 

< ay < Z 3 - - 1, j = 2,...,d, 

and periodic boundary conditions in the directions ej , 2 < j < d. Evidently (|2.47|) 
is derived from (|2.45|) . It is easy to see that the operator £ is invertible on the 
space L°°(fl) if the boundary conditions l|2.47() are imposed. In fact the solution to 
the equation 

(2.48) C V(x) = f{x), xetl, 

with boundary conditions l|2.47|l can be represented as an expectation for a continu- 
ous time Markov chain X(t), t > 0, on tt. For the chain the transition probabilities 
at a site x £ Vt satisfying < X\ < L\/2 — 1 are given by x — > x + ej, x — > x — ej, 
2 < j < d, each with probability l/2d, with periodic boundary conditions in di- 
rection e^, 2 < j < d. In the direction ei then x — > x + ei with probability 
l/2d + b(x) and x — > x — ei with probability l/2d — b(x). The waiting time at 
site x is exponential with parameter 1. If x\ = then x — > x ± ej, 2 < j < d, 
with probability l/2d[l + l/2d — < l/2d, and x — > a; + ei with probability 

[l/2d + 6(x)]/[l + l/2d — 6(a;)] < 1/d. The waiting time is exponential with pa- 
rameter [I + l/2d — b(x)]. Note that there is a positive probability that the walk 
will be killed at a site x with x\ = 0. A similar situation occurs at a site x with 
xi = Li/2-1. Nown z-ei with probability [l/2d-6(x)]/[l + l/2d+fe(x)] < 1/d 
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and the waiting time is exponential with parameter [1 + l/2d + b(x)]. The solution 
W of (|2.48(1 with boundary conditions (|2.47(l has the representation 



(2.49) = E 



f(X(t))dt\X(0) = x 



o 



x € 0, 



where t is the killing time for the chain. 

We may also consider the operator C of l|2.46[) with symmetric boundary condi- 
tions, 

(2.50) 

^(-l,x 2 ,...,x d ) = $>(0,x 2 ,...,x d ), y(L 1 /2,x 2 ,-,x d ) = %(Li/2-l,x 2 ,...,x d ), 

0< Xj <Lj-l, j = 2,...,d, 

corresponding to (|2.44|l . This is also associated with a continuous time Markov 
chain X(t) on Cl. The transition probabilities and waiting time at a site x £ Cl with 
< X\ < Li/2 — 1 are as for the chain defined in the previous paragraph. For x € Cl 
with x\ = reflecting boundary conditions corresponding to (|2.50(1 are imposed. 
Thus the waiting time at x is exponential with parameter [1 — l/2d + b(x)], x — > 
x ± ej, 2 < j < d, with probability l/2d[l — l/2rf + b(x)) and x — * x + ei with 
probability [l/2c? + 6(a;)]/[l — l/2d + b(x)]. A similar situation occurs atiefi with 
Xi = Li/2 — 1. The formal adjoint C* of the operator C of 1)2.46(1 is given by 

d 

(2.51) = VW-^ — Wx + ^ + Vix-ej)] 

3=1 

- b(x - ex)^(x - ei) + b(x + e x )*(a; + e^, ieii. 

It is easy to see that for functions <f>, on f2 satisfying the symmetric boundary 
conditions l|2.50|l there is the identity 

(2.52) ($£**> a = <*/:*) Aj 

where (-)a is the uniform probability measure on ft. Note that to show |2~52} one 
has to use the fact that the function b satisfies the antisymmetric conditions (|2.47|l . 
Hence the adjoint of the operator C acting on functions ^ : Q — > R with symmetric 
boundary conditions l|2.50|l is the operator C* of l|2.51|l also acting on functions 
with symmetric boundary conditions. In particular, it follows from Lemma 2.4 
that the solution if* of (|2.1t)|) . restricted to Q, is the unique invariant measure for 
the Markov chain X(t). 

Next let ipo '■ Cl — > R be the solution of the homogeneous equation (|2.48|l i.e. 
/ = 0, with the non-homogeneous antisymmetric boundary conditions 

(2.53) 

tpo{-l,X2, -,Xd) = -^0(0,^2, —,x d ), ip (Li/2,x 2 , -,x d ) = l-ip (L 1 /2-l,x 2 , ■ 

0< Xj <Lj-l, j = 2,...,d. 

One can see that i/jq is a positive function since it has a representation given by 
(|2.49fl . where / is the function 

f{x) = ^ + b(x),x£Cl, Xl =Lx/2-l, 
= 0, otherwise. 
The following lemma now shows that there is strict inequality in l|2.28|l 
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Lemma 2.6. Let ip* be the solution of 12. Iffy and ip be given by \2.11\j . Then there 
is the identity, 



(2.54) 



- + 2(^}=L?(^(.) 



55 



V'oO)xo(-) 



where xo '■ " ~~ > R is defined by Xo{%) = 1 if Xx = 0, Xo{ x ) = 0, otherwise. 

Proof. Since both tp* and ip are symmetric on ft in the sense of (|2.44H we may 
regard them as functions on ft with symmetric boundary conditions 12.50JI . We 
define a function tpi : ft — > R by rpi(x) = [Li/2 — 1/2 — xi]</j(x), a; € where tp 
is the solution to (|2.15() . ft is easy to see that 

(2.55) Ctpi(x) = [Li/2-l/2-xi]b(x)+tp(x), x G Cl, < x x < L x j2 - 1. 

We impose now symmetric boundary conditions on ■0i at x\ = 0, Xi = £i/2 — 1. 
One sees that (|2.55|) continues to hold at x\ — Li/2 — 1 but at xi = there is the 
formula, 

(2.56) AMa) = [ii/2 - l/2]6(x) + ^(x) - £i[l/2d- 6(x)]^(x). 

In deriving i|2.56[l we have used the fact that ip satisfies antisymmetric boundary 
conditions at X\ = 0. Now from (|2.16|) . Ij2.52jl . I|2.55jl . I|2.56j) we have that 



(2.57) 2(^V)=2(^»o = 



2L 1 {p*(-) 



5T*> 



<p(-)Xo(-) 



(<p*(-)[L 1 -l-2x 1 ]b(-)) i 



Next we define a function ip2 ■ " — » R by ip2(x) = (L\ — 1 — 2xi) 2 , x € fi. Then 
we have 



(2.58) 



A^fa) = 8[Li - 1 - 2xi]6(x) - 4/d, x e fi, < xi < Li/2 - 1. 



Again we impose symmetric boundary conditions on ip2 at xi = 0, x\ = L\j2 — 1, 
in which case ()2.58|l continues to hold at Xi = Zi/2 — 1. At Xi =0 there is the 
formula 

(2.59) Cip 2 {x) = 8[Li - 1 - 2x 1 ]b(x) - 4/d+ 4Li[l/2d- b(x)}. 

It follows now from l(2~TH| . lj2~5l>|) . J23HI, <|2~59l) that 



Xo(-) 



(2.60) -(^(•)[ii-l-2x 1 ]6(-))^ 
where we have used the fact that (<P*)q = 1. It follows now from i|2.57|) . (|2.59|) that 



(2.61) 



l/2d + 2{ip*ip) 



¥>*(•) 



1 

2d 



&(• 



[i+M-)]xo(-) 



We put now ipo(x) = [2xi + 1 + 4<^(x)]/2Zi, and it is easy to verify that ip 
satisfies the homogenous equation (|2.48() with the boundary conditions (|2.53() . The 
result follows then from 1|2.61[) . □ 

Proof of Theorem 1.1. The proof proceeds identically to the proof of Theorem 1.1 
of 6 , on using lemmas 2.1-2.4. □ 

Finally we wish to show that Theorem 1.2 holds to leading order in perturbation 
theory. 
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Theorem 2.1. There exists S > such that ifb : £1 — > R satisfies < sup w6 Q 1^(^)1 < 
<5 i/iera q(&) < l/2d) provided d = 1, or <i > 1 and Li < 4. If d > 2 and L\ > 6 f/ien 
i/iere exists arbitrarily small b with q(b) > l/2d. 

Proof. We shall use the LHS of H2.54JI as an expression for q(b). If b = then 
= 1, y> = =>■ ip = 0. Thus to obtain an expression for which is correct to 
second order in perturbation theory we need to expand ip* to first order in b and 
ip to second order. We consider first tp* which is the solution to 12.16fl . Letting 
A be the finite difference Laplacian acting on functions ^ : Q — > R with periodic 
boundary conditions, 

d 

A*(x) = ^2 + e i) + *( x - e i) _ 2 *(x)] , x G ft, 
i=i 

we have from <|2.5f II that 12. f 6fl is given by 

(2.62) -^ip*(x) + b(x + e 1 )i P *(x + e 1 )-b(x-e 1 )ip*{x-e 1 ) = 0, a; € 0, (99*) = 1. 

Since ( [r_ ei — r ei ] 6) = the solution to (|2.62l) is to first order in perturbation 
theory given by 

(2.63) if* = I + (-A/2C?)- 1 [ T _ ei - r ei ] 6. 
From (|2.46ll equation 12. f 5(1 is the same as 

(2.64) -—<p(x)-b(x)[T ei -T- ei ]p(x)=b(x), x e O. 
Using the fact that 

(6>-<6[r ei -r_ ei ](-A/2d)- 1 6) = 0, 
we see that the solution to (|2.64() correct to second order in b is given by 

(2.65) if = (-A/2d)~ 1 b + (-A/2d)~ 1 b [t 6i - T_ e J (-A/2d) _1 6. 

From l|2.f 7|l and Ij2.65|l we can obtain an expression for ijj which is correct to second 
order in b, 

(2.66) V = ^ [rei - r- ei ] (-A/2d)- 1 6 + 6 [r ei + r_ 6l ] (-A/2d)" 1 6 

+ 1 [r ei - r_ ei ] (-A/2d)- 1 & [r ei - r_ e J (-A/2d)" 1 6. 

From l|2.63|l . (|2.66|) we see that the lowest order term in the expansion of (ip*ip) in 
powers of b is second order. Thus correct to second order we have 

(2.67) = (b [r ei + r_ e J (-A/2rf)- 1 6) 

+ ^ (M^ - r_ ei ] (-A/2^)- 1 [r ei - r_ ex ] {-A/2d)~ l b) . 

The RHS of l|2.67[) is a translation invariant quadratic form, whence it has eigen- 
vectors exp[i£ • x], x G fl, with corresponding eigenvalue given by the formula, 
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We obtain an expression for the quadratic form 1)2. 67|) by doing an eigenvector 
decomposition in the x\ direction. Putting L = Li/2 we have that £1 = nk/L, 
k = 0, ±1, ±(L — 1), L. The function b then has a representation, 



If we use the antisymmetry property of b, b(y) = —b{2L — 1 — y) then one has that 

2L-1 L-l 

J2 b(y)e-^y = -2ie^' 2 £ b(y) sin&fo + 1/2). 

2/=0 y=0 

We conclude from this and 1)2. 68|) that the expression (|2.67|1 is the same as 



(2.69) r . 



Ad 
T 2 



E(-l)^(y)J + (E(-l)^(y) j ) 

V=0 I \y=0 ) I 

+^J2(( 1 Eb(y)^Hy + l)/L] 

k=l \ \y=0 / 

| cos ^ - 2 sin 2 ^ [-A d _i + 2(l - cos(7rfc/L))] ~ X j 

[-A d _ 1 + 2(l-cos(^/L))]" 1 ^6(y) S m7rfc(y + i)/Lj ^, 

where A^-i denotes the d — 1 dimensional Laplacian acting on the space {x\ = 0}. 

Observe now that the L dimensional vectors sin7rfc(y + l/2)/L,0 < y < L — 1, 
are mutually orthogonal, k = 1, ...,L. This is a consequence of the fact that they 
are the eigenvectors of the second difference operator on the set {0 < y < L — 1} 
with antisymmetric boundary conditions. It follows that the quadratic form (|2.69() 
is negative definite if and only if all the eigenvalues (|2.68l) are negative. This is the 
case for d = 1. For d > 1 it is still true provided L < 2, but already for L = 3 it 
is false. Thus for L — 3 one can find a b such that the homogenized limit has an 
effective diffusion constant which is larger than the 6 = case. □ 

3. Proof of Theorem 1.2 

We shall use the representation for the effective diffusion constant given by the 
RHS of l|2.54|l . We consider first the d = 1 case. 

Lemma 3.1. Let O be the space fl — {xeZ: l<x< L}. If ip* : — > R is the 
solution to $2.16]) then y*(l) is given by the formula, 

L L r-l L 

(3.1) <p'(m=L nvE if* n ^ 

k— 1 r— 1 j=l j"=r+l 

where the 5j, 5j 7 1 < j < L 7 are given by 

(3.2) 6 j = 1/2 -b(j), S j = 1/2 + b(j). 
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Proof. From (|2.51(l we see that f* : f2 — * R satisfies the equation 

(3.3) cp*(x) - - tp*(x + 1) + - 1) - 6(o; - l)<p*(x - 1) 

+ &(z + l)v?*(x + l) = 0, l<x<L, 
with the symmetric boundary conditions and normalization given by 

(3.4) p*(0) = V *{L + 1) = ip*(L), = 1. 

We can solve (|3.3f) . (|3.4|) uniquely by standard methods. Thus putting Dip*(x) = 
<p*(x) — ip* [x — 1), 1 < x < L, then we may write l|3.3|l as 

(3.5) i Dtp*(x)-Dtp*(x+1) -b(x-l)<p*(x-l)+b(x+l)<p*(x+l) = 0, 1 < x < L. 
If we sum l|3.5|l over the set {I < x < y} we obtain the equation, 



-[2V(1) - JV(l/ + 1)J - &(0)¥>*(0) - 6(1)^(1) 

+ &(y)p* (y) + % + 1 V*(y + 1) = 0, I < y < L. 
Then, using the fact that — ¥?*(0), b(l) = —6(0), we conclude that 

<p*(y + l) = S y <p*(y)/5 y+1 , l<y<L, 

whence we have 

(3.6) <p*(y) = <^(1) J] A,-/5 i+x , I < y < L. 

3=1 

The formula l|3.1|l follows from (|3.6() and the normalization condition in (|3.4() . □ 

Lemma 3.2. Lei ft be the space ft — {x € Z : 1 < x < L}. If i))q : Cl — ► R is the 
solution to the homogeneous equation {2.4tf) with the boundary conditions \2.53\) 
then ipo (1) is given by the formula 

L L r-l L 

(3.7) 2^o(i)=n 5 VE ik n 

fe=l r=l j'=l j=r+l 
where Sj, Sj, 1 < j < L, are as in $3.2\) . 

Proof. From (|2.46l l|2.48() . I|2.53|l . we see that ipo(x) satisfies the equation, 



(3.8) M^)-^[Mx+l)+Mx-l)\ -b(x)\ip (x+l)-Mx-±) 



0, I < x < L, 



with the boundary conditions, 

(3.9) ^o(O) = -Vo(l), ML + 1) = 1~ML). 

We can solve (|3.8f) . (|3.9|l by standard methods. Thus putting Dipo(x) = ipo(x) 
i/)o(x — 1) equation l|3.8|l implies 

(3.10) Dip (x + l)=S x DilJo(x)/S x , l<x<L. 
Observing from l|3.9|l that Dipo(l) = 2ip (l) we see from (I3.10|l that 



(3.11) 



DMv + 1) = 2VVi) II 5 i/ s J> 1 <y< L - 

3=1 
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If we sum (|3 . 1 1|) we obtain a formula for ipo(y) given by 

{y r-l 1 
r=2j=l J 

Since Dip (L + 1) = 1 - 2ip Q (L) from (|3~3|) it follows that if we add (f5~TT|) to twice 
(|3.12[1 when y = L we obtain a formula for V'o(l) given by 

L ( L L r-l L L 

(3.i3) 2Mi) = n v \ n + 2 e n ^ n + n ^ 

fe=l t- 7 ^ 1 r=2j=l j=r j=l 

One can easily see that the denominator of the expression in (|3.13|) can be rewritten 
as in (E7J). □ 

Remark 1. Observe from that under the reflection b — ► —6 the expres- 

sion tp*(l)5i/L becomes 2ipo(l). 

Lemma 3.3. There is the inequality <^*(l)<5i-0o(l) < 1/8L. 
Proof. For 1 < r < L define a r by 



r-l L 



if* n ^ 

j=l j=r+l 



and a r the corresponding value of a r under the reflection b — > —b. From Lemma 
3.1, 3.2 we see that we need to prove that 



4L 2 ]JS k 6 k < \j2 a r] fe^f- 

k=l [r=l ) [r=l J 



Using the fact that for 1 < r, s < L, 

(a r a s + a s a r )/2 > (arOrasCis) 1 / 2 , 

we see that 



a r > < a r > > (ararasds) 1 / 



1/2 

% r— 1 J { r—1 ) r,s—l 

L L L 



r,s=l k=l k=l 

where we have used the fact that for 1 < j < L, one has 5jSj < 1/4. □ 

Proof of Theorem 1.2. {d = 1): This follows from Lemmas 3.1 - 3.3 and Lemma 
2.5, using the RHS of (12.54(1 as the representation for q(b). □ 

Next we turn to the d > 1 case with L\ — 2. Then we can write Ct = {(0, y) : 
y € f^d-i} where £ld-i C Z is a d — 1 dimensional rectangle. It is easy to see 
now from (|2.ti2|) . on using the anti-symmetry of b and the symmetry of ip*, that 
ip* = 1. Also from (|2.64|) . on using the anti-symmetry of b and tp, we have that 



(3.14) ^(0,y) = 2d[-A d _ 1+ 4] ^(O.y), 
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where in Q3.14|l the operator Ad-i is the discrete Laplacian acting on functions 
* : Qd-i -> R- We have then from f2~TT)l that ip(0,y) = -2b(0,y)(p(0,y), and so 
we get the formula for the effective diffusion constant, 

(3.15) l/2d+ 2 < cf>*il> >= l/2d-8d(&(-) [-A d _ x + 4] _1 &) 

It is clear that the RHS of l|3.15[l is smaller than l/2d. We can alternatively derive 
the effective diffusion constant formula by using the expression on the RHS of H2.54JI . 
Thus we have 

MO, V) = 2^ [-Ad_i + 4]- 1 [l/2d + 6(0, y)], 
whence the effective diffusion constant is given by the formula 

(3.16) Ll{p(-)[l/2d-b(-)]M-)Xo(-)h = 

8d([l/2d -&(•)] [-A d _! +4]- 1 [l/2d + &(•)]) 

We shall use the formula on the LHS of (|3.16(l to obtain an expression for the 
effective diffusion constant in the case L\ = 4. Here il is the space (l = {(n,y) : 
n = 0, 1, y G f2(j_x}. For y € fi<z-i we define S y , S y , e y , e y by 

(3.17) 6 y = l/2d-b(0,y), S y = l/2d + b(0,y), 
e y = l/2d + 6(l,y), e y = l/2d-b(l,y). 

We see then from (|2.fi2|) . H3.17(l that <p* satisfies the system of equations, 
(3-18) ( ~ Ad 2 d +2 ) y * (0 ' y) " e y tp*(l,y)-5 y( p*(0,y)=0, 

i ~ Ad 2d +2 ) v * {hv) " £ ^*(i,y)-V*(o,y) = o. 

Adding the 2 equations above we conclude that — Ad_i[y*(0, y) + y>*(l,y)] = 0, 
y G f^d-i, whence on using the normalization < ip* >^= 1 we conclude that 
<p*(0,y) + <p*(l,y) = 2, yGfirf-i. Hence from (|3.18|l we have that y*(0, ^ satisfies 
the equation, 

(3.19) [-A d _i/2d + e y + S y ] V *{0,y)=2£ v , y G ft d _i. 

Evidently IpTTOjl has a unique positive solution. 

We proceed similarly to obtain a formula for ipo . Thus -0o satisfies the system of 
equations, 

(3.20) ( ~ Ad " 1 + 2 Wo,y) - ^^ o (l,y) + ^o(0,y)=0, 



2d 

-A d _! + 2 



■00(1,2/) + £ y 1po0-,y) ~ £y^o{0,y) = £y, 2/ € fid-1. 



2d 

Adding the two equations in l|3.20|l we get 

(3.21) [-A d _ x /2d+e y + 6 y ] {V»o(0,y)+Vto(l, !/)}=£», V e fid-i- 
We may also rewrite the first equation of l|3.20|l as, 

(3.22) ( Ad 2 d + 4 ) M°>v) = ~ 5 viM0,y) + Mhy)}, i/eo^ 
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We conclude from ijOty . lj3~2T)l . and the formula on the LHS of (|3~ll)|) that 

the effective diffusion constant q(b) is given by, 

(3.23) q(b) = 2 7 d 3 ( [d [-Ad-i + 2 - 1/p 1 e} 

[-A d _i + 4]' 1 f I [-A d _i + 2 + V]" 1 e} ) 

^ J I fid— 1 

where V : fid-! ^ R is given by V(y) = 2d[&(l,y) - 6(0, y)], y € 
We first show that q(b) < l/2d in the case where V is constant. 

Lemma 3.4. Let q(b) be given by \S. 23(1 and assume V is constant. Then there is 
the inequality, q(b) < l/2d. 

Proof. Since e + S = (2 + V)/2d there is an / : ttd-i — > R such that 

(3.24) e = (2 + V)/4d+/, * = (2 + V)/4d - f, 
e = (2-V0/4d-/, J=(2-F)/4d+/. 

We rewrite the expression on the RHS of l|3.23|l in terms of /. To do this we let 
w+,W- be solutions to the equations, 

(3.25) [-A d _ 1 + 2 + V)w+ = f, 

[-A d _ x + 2 - V] w- = /. 

It follows that 

[-A d _i + 2 + VY l e = l/4d + w+, 
[-Ad-i + 2-Vr 1 e = l/M-w-. 
Hence from (j3~2Hl q(b) is given by the expression, 



q(b) = 2 7 d 3 ' 



V 



Ad 



f 



Ad 



[-A d _ x + 4]- 



V 



Ad 



+ f 



1 

Ad 



W-i 



This is a quartic expression in / and the zeroth order term is given by, 
(3.26) zeroth order = ^- ((2 + V) [-A d _j + 4] _1 (2 - V)) 

Observe that the expression in (|3.26|) is identical to the RHS of l|3.16|l if e = S. For 
the first order term we have the expression, 

(3.27) 

first order = 2 ((2 + V) [-A d _ x + 4]" 1 /) - 2 (/ [-A«j_i + 4]" 1 (2 - V)) 
+2 ((2 + V) [-A d _! + 4]" 1 (2 - V>+) n "2 ((2 + V>_ [-A«,_i + 4]" 1 (2 - 7)) 
Now from l|3.25[) we have that 

(3.28) (2-V>+ = [-A d _i + 4] w + - /, 

(2 + V)u>_ = [-A d _i + 4] w_ - /. 



HOMOGENIZATION OF A DIFFUSION PERTURBED BY A VECTOR FIELD 

From this we conclude that the expression in Ij3.27|l is the same as 

2 ((2 + V)w+) ttd _ i - 2 <(2 - V)w-) Qd _ i = 0. 
The second order term in l|3.23[l is given by 
(3.29) 

-2 3 d If [-A d _x + 4]- 1 /) -2 3 d ((2 + V>_ [— A<j_i + 4]" 1 (2 - V)w, 
- 2 3 d (/ [-A d _! + 4]- 1 (2 - V>+\ - 2 3 d ((2 + V>_ [-Aj_i + 4]" 1 f) 

+2 3 d ((2 + V) [-Ad_i + 4]" 1 fw+) +2 3 d (fw- [-A d _! + 4]" 1 (2 - F)\ 

Observe from 1|3.28[1 that there is the identity, 
(3.30) 

((2 + V>_ [-Aj_i + 4]- 1 (2 - V>+) = (/ [-A d _ a + 4]" 1 /) 

We similarly have that 
(3.31) 



/ [— A<j_i + 4p (2 - V)w+) = - (f [— Ad_i + 4p /) + (/«;+>„ 

((2 + V>_ [— Ad_i + 4]- 1 /) = - (/ [-A rf _ x + 4]- 1 /) + (/„,_> 

Define now [7 : £7^—1 ~> R- as the solution to the equation, 

(3.32) [-A d _i + 4] U = V. 

Then from equations H3.29|l - (|3.32(l we see that the second order term in (|3.23|l 
given by 

(3.33) second order = 2 3 dj - (w- [-A d _ a + 4] w + ) Qd _ i 

+ \ (/ [w + + + (fU [w + w_]) nd i | . 

The term of third order is given by 
(3.34) 

2 5 d 2 If W - [-A d _! + 4]- 1 /) +2 5 d 2 (/ W _ [-A d _! + 4P 1 (2 - V)w+) 

- 2 5 d 2 (/ [-Ad_i + C 1 fw + )- 2 5 d 2 ((2 + V>_ [-A rf _! + 4]" 1 . 

Using l|3.28|l again we see that the expression (|3.34|) is the same as 

2 5 d 2 </ W - w+) nd i - 2 5 d 2 </ w- w + ) nd i = 0. 
Finally the fourth order term is given by 
(3.35) fourth order = 2 7 d 3 (/ w_ [-A d _i + 4] _1 / w A 



Hence from (|3.26ll . (|3.33|l . (|3.35() we have the formula for q(b), 
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(3.36) 

q(b) = ^{( 2 + V ) i- A d-i + (2 - V))^ i - 2 3 d («,_ [-A d _ a + 4] W+ ) n<j _ i 
+2 2 d(/[ W+ + W _])^ i +2 3 d(/C/K- W _]) 0d i +2 7 d 3 (/ W _[-A d _ 1 +4]- 1 /w+) 

\ / Od-i 

It is clear from the definitions of V, / that 

(3.37) \V(y)\<2, \f(y)\<[2-\V(y)\]/4d, y G fW 
From (|3.37|) it follows that there is the inequality, 

(3.38) 2 7 d 3 (/ W _[-A d _ 1 +4]- 1 /^+) <2 4 d 3 ((/ w _) 2 ) f2 

\ l (ld-1 '"d-l 

+2 4 d 3 ((Z^) 2 )^ < d([2 - \V\} 2 [w 2 _ + wl})^ . 
We define now a quadratic form Qy depending on V by 



(3.39) Q v (f) = (w- [-A d _ x + 4] w + ) nd i - - (f[w + + w-])^ 

- (fU[w + w.}) Qd _ t - I ([2 - \V\] 2 [ W 2 _ + wl]) nd i . 
It is evident from (|3~3(j|) . [|3~5gj> that 

(3.40) q(b) < i ((2 + K) [-A d _x + 4]- 1 (2 - 10) i - 2 3 d Qy(/). 

Thus to prove the result it will be sufficient to show that Qy is nonnegative definite. 

Since V is constant we can compute the eigenvalues of Qy explicitly. Thus if p 2 
denotes the eigenvalue of — A^-i, corresponding to the eigenvector exp[i£ ■ x],x £ 
fid-i, then the eigenvalue of Qy is 

(3.41) 1^2(4 + V 2 )(p 2 + 2 + V){p 2 + 2-V)- 

i(2 - \V\) 2 [{p 2 + 2 + V) 2 + {p 2 + 2- Vf] | J 4(p 2 + 2 + Vf{p 2 + 2- V) 2 . 
We can rewrite the numerator of (|3.41|) as 

(3.42) {2(4 + V 2 ) - (2 - l^l) 2 } [p 2 + 2] 2 - V 2 {2(4 + V 2 ) + (2 - \V\) 2 } . 

Since \V \ < 2 the expression in H3.42|l is bounded below by its value for p = which 
can be written as 

{4 + V 2 ){2- |V|)[2 + 3|V|] > 0. 

□ 

We proceed now to the general case which will follow from the fact that the 
quadratic form l|3.39|) is nonnegative definite for any V satisfying (|3.37|l . From here 
on we shall denote A<2_i,(-) n i simply as A and (•) respectively. We first note 
that by using (|3.25|) we can obtain some alternate representations for Qy. Thus if 
we write 

(ui_ [— Ad_i + 4] w+) = 2 (w- w+) + - {f[w+ + w-]) , 
we see that Qy is given by 

(3.43) Q v (f) = 2 {w- w + ) - (fU[w + - «;_]) - g ([2 - \V\] 2 [w 2 _ + w\]) . 
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We also have that 

(fU[w+ - w-}) = ({U[—A + 2- V]w-} w+) - ({C/[-A + 2 + V]w+} W-) 

= -2 (UVw- w+) - (Uw + Aw-) + (Uw- Aw+) . 

Hence from ()3.43|) we have the formula, 
(3.44) 

Qv(f) = 2 (w- w+[l + UV]) + (Uw+ Aw-)-{Uw- Aw + )-~ <[2 - \V\] 2 [w 2 _ + w\]) . 

8 

We first show that a simple quadratic form related to Qv is nonnegative definite. 

Lemma 3.5. Let V satisfy jff.^TI ) and w + , w_ be solutions to {K 25\) for any 

f : £ld-i ~~ * R- Then there is the inequality (w+ W_) > 0. 

Proof. Let w be the solution to the equation, 

(3.45) [- A + 2 + V] [-A + 2- V]w = f. 

Then from (|3.25|) we see that w+ = [—A + 2 — V]w. Hence we have from (|3.25H . 
(l3~4*5jl that 

(w+ w-) = ([-A + 2 - V]w w-) = (w[-A + 2 - V]w-) 

= (w f) = ({[-A + 2 + V]w}{[-A + 2-V] w}) 

= ({{-A + 2)w} 2 ) - (V 2 w 2 ) > ((4 - V 2 )w 2 ) > 0. 

□ 

To proceed further we need to localize the quadratic form l|3.44ll . 

Lemma 3.6. Let £ + , £_ be operators on functions $ : Qd-i — * R defined by 

£+$ = (-A + 2)$/V-$, 
£_$ = (-A + 2)$/V + $, 

where we assume V satisfies \3.3 r ]\l and V(y) / 0,;/ € Qd-i- Then C+,C~ are 
invertible and there is the identity, 

(3.46) [-A + 2 + V]£+ = [-A + 2- V}£- . 

Proof. Verification. □ 

It follows from Lemma 3.6 that we can choose / in l|3.25[l . (|3.44|) as the operator 
(|3.46Jl acting on a function $ : ild-i —> R, in which case w + = C+Q, w- = 
If we substitute into (|3.44|) we obtain a quadratic form Qy($) which is local in 

and it is this quadratic form which we will show is nonnegative definite. First 
we show that the quadratic form obtained from Qv upon replacing U by V/4 is 
nonnegative definite. 

Lemma 3.7. For $ : fld-i — > R and w + = w- = there is for d = 2 

the inequality, 

(3.47) 2 (w_ w+[4 + V 2 ]) + (Vw+ Aw_) - (Vw- Aw + ) 

^([2-\V\} 2 [w 2 _+w 2 + })>0. 
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Proof. We first note that the first term in (|3.47|) is nonnegative. Thus we have 

(w. w+[A + V 2 ]) = ([(-A + 2)$] 2 [4/V 2 + 1]) - ($ 2 [4 + V 2 ]) 

> 2 ([(-A + 2)$] 2 ) - ($ 2 [4 + V 2 ]) > ($ 2 [4 - V 2 ]) > 0, 

where we have used (|3.37|) . The second and third terms of l|3.47[l are given by the 
formula, 



(-A + 2)$ 



V 



(Vw + Aw_) - (Vw_ Aw + ) = 

2 ([(-A + 2)*](A $)) - 2 ([A(V$) 
on summation by parts. From the last two equations we therefore have that 
(3.48) 2 (w- w+[A + V 2 ]) + (Vw+ Aw_) - (Vw- Aw+) 



8 



[(-A + 2)$] 



V 2 

Using the fact that 
1 



2(V 2 <f> 2 ) +4((V$) 2 ) -2([A(V$)] 



(-A + 2)*] 



V 



V(x) 



A(V$)(x) = 2(d-l)$(x) 



d r 

£ 

J=2 



V(x + ej) V(x-ej) 



V(x) 



V(x) 



$(x - ej) 



we conclude from l|3.48|l that 



(3.49) (w- w+[4 + V 2 ]) + - (Vw+ Aw-) - - Aio + ) 



[(-A + 2)d>] 
V 2 



-\ + ([4(d - 1) - U 2 ]$ 2 ) + 2d <(V$) 2 ) 



[(-a+2)$(.)]x; 

» J=2 



-*(-+e,- 



Now the Schwarz inequality yields 
V(»)*(lO 



[(-A + 2)$(a;)] 



y(x) 



V(-) v J ' V(-) 

$(y) 2 1 [(-A + 2)$(x)] 2 



V(y) 2 Aa V(x) 2 

x,y € ftd-i, 



^ _ . 2 1 [(-A + 2)$(x)] 2 
< a$(y) 2 + 1 J 



for any a > 0. Hence there is from 



1 



a V(x) 2 
the inequality, 
1 



(3.50) (w- w+[4 + V 2 ]) + - (Vw+ Aw_) - - (Vw- Aw+) 



> 



2(d- 1) 



-A + 2)3f 
V 2 



( [2(d - 1)(2 - a) - U 2 ] $ 2 ) + 2d <(V$) 2 ) 
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For d = 3 and a = 1 the RHS of H3.50|) is evidently nonnegative but this is no 
longer the case when d > 3. For a = 1, d = 2 the RHS of H3.50JI is bounded below 
by the nonnegative quantity, 

(3.51) ~ ( (^^) [(-A + 2)$] 2 ^ + <(4 - T/ 2 )$ 2 ) . 

We consider now the last term in the expression (|3.47|l . We have that 

(3.52) \ <[2 - |U|]> 2 + = \([2- \V\? { [( ' A + 2)$]2 + <t> 2 

If we now use the inequality [2 — IV^] 2 < 4 — V 2 we see that the expression (J3.52JI 
is less than (13.51(1 . Hence the inequality l|3.47jl is established for d = 2. □ 

Next we turn to showing that Qy is nonnegative definite for d = 2. To do this 
we write the solution U of l|3.32|l as 

(3.53) u( x ) = ^j2 G ^ v ( x + y^ 

y 

where G{y) is the Green's function for [— A/4 + l] -1 , whence G(y) is nonnegative 
for all y and 

(3.54) ^G(y) = l, G(y) = G(-y), y = 1,2,... . 

We consider the first three terms in the expression H3.44|l for Qv- Using Lemma 
3.6 and setting w+ = = we have that 

(3.55) (w-w+[l + UV}) + i (Uw+Aw-) - i (Uw-Aw+) 
[(-A + 2)$] 2 [1 + £7V] ) - ($ 2 [1 + UV]) 



V 2 

^[A$] [(-A + 2)$A - /A(C/$)i[-A + 2]$ 
[(-A + 2)$] 2 ^-($ 2 )-(<l>W) 
2 ^*(-A + 2)$^ - ^A([/$)i[-A + 2]$ 

^[(-A + 2)$] 2 ^-($ 2 )-($W) 
4 /~*(-A + 2)$^ - /( ri f/)( Tl $)I[-A + 2]$ 

(r_il/)(r_i*)i[-A + 2]$V 

where r x tp(y) = (p(x + y), yeZ. We consider the last three terms in the previous 
expression. We write using H3.53fl . 

4 (£$(-A + 2)$\ = G(0) <(V$) 2 + 24> 2 ) + ^ G(y) *(-A + 2)*\ , 
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■I>)-i-A-2>\ = ~G(-1){(V(ti$)V$> + 2 <(>!$)$)} 



V 



I^G( y -l)(^Ir 1 $(-A + 2)$ 



with a similar expression for the last term in (|3.55|) . We conclude from this that 
the last three terms of (|3.55|) are given by, 

(3.56) 

4 ( ^$(-A + 2)$\-/( n C/)(n$)^[-A + 2]#\-/(t_iC0(t_i$)^[-A + 2]$ 



1 



1 



G(0) <(V$) 2 + 2$ 2 )--G(-l) {(V(t!$)V$> + 2 ((n$)$)}--G(l) {(V(r_!$)V$) + 2 ((t_i*)$)} 



+ E 
E 

8,<-l 



G(y) - \ G(y - 1) - \ G(y + 1) 
G(y) - ~ G(y - 1) - \ G(y + 1) 



^- n$(-A + 2)$ 

^ T_!$(-A + 2)$ 



E G ^)(lT [*-n$](-A + 2)$)-i^G(y+l)(^ [r_ 1 <i>-T 1 $](-A + 2)<i> 



1 

4 



V 



E G ^)(^ [*-r- l $](-A + 2)$)-\ Yl G(y-l)(^f [r 1 d>-r_ 1 $](-A + 2)$ 



r 



y<-l * ' y<-l 

We shall use the representation (|3.56(l to show that the quadratic form Qy is non- 
negative definite. 

Lemma 3.8. Suppose the function G(y) of if ff. 53\) is decreasing, non-negative for 
y > 1, satisfies \3. 54{ l and the inequalities, 

(3.57) (-A + 2)G(y) < 0, y > 1; 1 - G(0) - 2G(1) < G(l)/2; G(2) < G(l)/5. 

Then the quadratic form Qv of \'d-4-4\l * s nonnegative definite. 

Proof. We estimate the terms in (|3.56[l by applying the Schwarz inequality. Before 
doing this we make one further simplification of terms in l|3.56[l . We write 
(3.58) 

E G (y)(^ [*-ri$](-A + 2)$\+^ G{y)hf [$ - T _ 1$ ] (-A + 2)$ 
y>i ^ ' v<-i N 

= E G ^)( I |r [*-ri$](-A + 2)*\+5^ G{3,)(^- [$-r_!$](-A + 2)$ 



y>2 

+G(i) ( n v 



r 



l v 

~ 7 



[$ - n$] (-A + 2)$ }+G(-l) ( r_iV 



G(1) ((tiV)V[* - n$]A$)-^-^ ((r_iV)^[4 - r_x$]A$) 



i y 

y ~ 7 

G(l) 



4 * ^ 4 sv ' 1 2 

where we have used the fact that G(l) = G(— 1). We also rewrite the first two 
terms on the RHS of H3.58J) as 

'(TyV-V) 



[$-t_i$](-A + 2)$ 

((^)nvi-) 2 ), 



(3.59) 



E G (y) 

y>2 



r 



[$-n$](-A + 2)$ 
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E G (y) 

y<-2 



(TyV - V) 



V 



[$ - t_!$] (-A + 2)$)+E G(v) <( A$ ) 2 + 2(V$) 2 ) . 
We similarly rewrite the sum of the last and third last terms of (|3.56|) as 



(3.60) 



y>l 



V 



[r_!$-ri$] (-A + 2)$ 



y<-l 



V 



[ n $ _ y^^] (_A + 2)$ 



Consider now the first three terms on the RHS of l|3.55[l . These can be written as 
(3-61) ( (?2 - l) [(~ A + 2 ^ 2 ) + I (( A$ ) 2 ) + ( (V$ ) 2 ) - 

i G(o) <v 2 $ 2 ) + i E G (y) [<W vf&) - ({{r y Vf + V*} $ 2 )] . 



Next we apply the Schwarz inequality to terms in H3.56J1 . Thus we estimate 



(3.62) 



^ Tl $(-A + 2)$ 



^ ($ 2 ) + (-^[(-A + 2)<i>] : 



with a similar estimate when ti$ is replace by r_i$. 
Observe now that 



E 

y>l 



G{y)~- A G(y-l)-- A G{y + l) 



1 



[2G(0) - 1 - G(l)], 



where we have used (|3.54f) . Hence on using the fact that [—A + 2]G(y) < 0, y > 1, 
we see from <|3.62ll that the sum of the first five terms on the RHS of 1|3.56[) are 
bounded below by the expression, 

(3.63) [G(0)-Jg(1)]((V$) 2 + 2$ 2 ) 



<* 2 > 



1 .[(-A + 2)$] 2 \l[G(0)-i-iG(l)]. 



2 >[(-A + 2 )$ ] 2 /J „ . . 2 2 

If we combine the estimate (I3.63|) with (|3.61|l and use the fact that | | < 2 we get 
a lower bound for the sum of the first three terms of (|3.55|l and the first five terms 
of (|3.56|l . It is given by, 

^-\) [(-A + 2)*] 9 ){§^G(1)-G(0) 



(3.64) 



+ \ ((A$) 2 ) |^ + 1 G(l) - G(0)| + \ ((V$) 2 ) + \ G(0) ((4 - V 2 )&) 
+ \ E ^ ~ ^) 2 $ 2 ) + [l-G(0)] < $ 2 > -\ E G(y) {{(r y Vf + V 2 } $ 2 ) . 

y^O y^O 

Observe that all terms in (|3.64|) except for the final one are nonnegative. Further- 
more, the sum of the last two terms is nonnegative. 

Next we estimate the terms on the RHS of <|3.58ll which involve G(l) and G(— 1). 
To do this we use the Schwarz inequalities, 

' { 1 [$-ti$] (-A + 2)$ 



(3.65) 



71 V 



V 
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a-i 



1 

V 2 



[(-A + 2)$] 



-<(V$) 2 ), 

Q'l 



\{{nV)V [$ - n$] A$) | < a 2 <(A$) 2 ) + — <(V$) 2 ) 

«2 



for any constants <X\,ol% > 0. Hence on using the fact that \V\ < 2 we see that the 
expression (|3.64|) plus the terms in G(l), G(— 1) of 13.58fl is bounded below by the 
expression, 

(3.66) 

1 



1 



1 



_-_)[(-A+2)$r 



| + ~ G(l) - G(0) - 2a x G(l) 



<(A$) 2 ) 



3 

2 ' 2 



^G(l)-G(0)-4a 2 G(l)}+((V$) 2 ){ 
- G(0) ((4 - y 2 )$ 2 ) + i £ G (J/) ((TyV - v?&) . 



|-2G(1)--G(1)-- 
2 ai a2 



Let us assume for the moment that G(y) = for y > 2. Then (|3.55|l is bounded 
below by J3~m We write G(0) = 1 -7 whence G(l) = 7/2. Since (-A + 2)G(y) < 
0, y > 1, we must have 7 < 1/3. We choose «i such that 



\ + \ G(l) - G(0) - 2aiG(l) = \, 



(3.67) 

which yields ai = 5/4. We choose a 2 so that 
(3.68) 



| + ^G(l)-G(0)-4a 2 G(l)=0, 



which yields a 2 = 5/8 + I/47. The coefficient of < (V$) 2 > in Ij3.66jl is therefore 
bounded below by 1.5 — 2.67 > since 7 < 1/3. Hence from <|3.52[1 the quadratic 
form Qy/2 of (|3.44|) is bounded below by twice the expression, 



(3.69) 1 _ ^ [(-A + 2)$] 2 ^ + \ G(0) ((4 - V 2 )<5> 2 ) 

[(-A + 2)$] 2 



-o <[2-|VT 



V 2 



If we now use the fact that [2 — | V |] 2 < 4 — V 2 we see that <|3.69[1 is nonnegative. 

To complete the proof of the lemma we need to estimate the sum of the terms 
in (|3~59"|) . IE2I0J). We rewrite these as 



(3.70) 



- 4 G(2) 



(nV-V) 



v 

{t-xV-V) 



[t_i$-ti$](-A + 2)$ 



V 



G(y)--G(y + 1) 



E 

y>2 



y>2 



(TyV - V) 

V 



E 

y<-2 



1 



G(y) --G(y- 1) 



[n$-r_i$] (-A + 2)$ 

(T " V ^ [ $ - r i $ K-A + 2)$ 
$-r_i$](-A + 2)$ 

[$-r_i$] (-A + 2)$ 
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y<-2 



(r y V - V) 
V 



[$ -ri$] (-A + 2)$ 



+ ^[1-G(0)-2G(1)]((A^) 2 + 2(V^) 2 ). 
We first estimate the third term in H3.70JI . Thus we write 

'(TyV-V) 



(3.71) 



V 



(Ty V-V) 



[$-ri$] (-A + 2)$ 



1 

V 



V 

1 



[$-n$](-A + 2)$ 



-~ ((r a V - V)V[$ - n$]A$) + ~ ((r, 1/ - V)V[$ - Ti$]$) . 

We estimate the first two terms on the RHS of l|3.71[l similarly to <|3.65ll . For 
the third term we use 

(3.72) | {(r y V - V)V[* ~ ri*]$) | < a ((V$) 2 ) + i (( Ty V - V) 2 ^ 2 ) , 

for any a > 0. Choosing a — 4 in (|3.72l) it follows that the sum of the third, fourth, 
fifth and sixth terms of (|3.70l) is bounded below by 



(3.73) 



|y|>2 

a 4 ((A$) 2 ) 



<(V$) 2 ) + i((r^-^) 2 $ 2 )|. 



2 1 
1 1 

for any 03, 0:4 > 0. We estimate the sum of the first two terms in (|3.7U|I from below 
similarly. Choosing now a — 2G(2)/G(1) in l|3.72|) we obtain the lower bound, 



(3.74) 



I 1 G(2) 



G(2) 

a 6 ((A$) 2 ) 



a 5 a 6 G(l) 



((V$) 2 ) 



- ^((nV-Vf^-^iiT^V-V) 2 * 2 ), 

for any CH5, ag > 0. 

We may now obtain a lower bound for (|3.55(l by adding l|3.66(l to the final term 
in (|3.7U|) and the expressions of (|3.73|) and l|3.74|l . We obtain the lower bound, 



(3.75) 



yS-i) [("A + 2)$] 



G(0) 



2aiG(l) - 2a 5 G(2) - 2a 3 [l - G(0) - 2G(1)] 



+ - ((A$) 2 ) - + - G(l) - G(0) - 4a 2 G(l) - 4a 6 G(2) 



4a 4 [l - G(0) - 2G(1)] + 2[1 - G(0) - 2G(1)] 



+ <(V$) 2 ) 



2G(1) 



ai 



-G(l) 



1 

a 2 



G(l)-G(2) 2 /G(l) 
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2 1 

a 5 a 6 



G(2)- 



2 1 
— + — + 2 

Ct3 Ct4 



[1 - G(0) - 2G(1)] 



+ [1 - G(0) - 2G(1)]\ + i<?(0) <(4 - U 2 )$ 2 ) . 
We may rewrite the coefficient of the first term in (|3.75[) as 



(3-76) - 
1 



- G(l) - G(0) - 2 ai G(l) - 2a 5 G(2) - 2a 3 [l - G(0) - 2G(1)] 



= - + (l-2a 3 )[l-G(0)-2G(l)] 



5 „ „ G(2) 
--2 ai -2a, W) 



G(l). 



If we set now 
(3.77) 

we see as in (|3.67|) that the coefficient of the first term in (|3.75[1 is 1/2. We similarly 
rewrite the coefficient of the second term as 



a 3 = l/2, ai + a s G(2)/G(l) = 5/4, 



(3.78) 



| + i G(l) - G(0) - 4a 2 G(l) - 4« 6 G(2) 
-4a 4 [l - G(0) - 2G(1)] + 2[1 - G(0) - 2G(1)] 



1 



= - + (3-4a 4 )[l-G(0)-2G(l)] + 



5 A A G ( 2 ) 

4a 2 - 4a 6 — ^4 

2 2 6 G(1) 



G(l). 



Hence if we set 

(3.79) a 2 = 5/8, a 4 = 3/4, a 6 = l/8G(2) 

then the second term in l|3.75|l is zero. We consider the third term in (|3.75(l . This 



can be written as 
3 



(3.80) 



_ 2G(1) - -G(l) - -G(l) - G(2) 2 /G(l) 
I a\ a 2 



' 2 


1 


G(2)- 


to 




1 


— + — + 1 


_a 5 


a 6 


a 3 a 4 



--(-A + 2)G(1) 



7-2—1' 

ol\ a 2 



2 1 

a 3 a 4 

2 



G(1)-G(2)7G(1) 



[1-G(0)-2G(1)] 
[1-G(0)-2G(1)] 

G(2). 



3 2 1 

2 as as 

Using the inequalities (|3.57|) we see this is bounded below by the expression, 



(3.81) 



7- 



2 


1 


1 


'1 


2 


1 " 


1 


1 


'3 


2 


1 

H 

a 6 


ai 


a 2 


h 2 


2 


a 3 


a 4 _ 


25 


5 




«5 



G(l) 



6.91 



G(l). 



2 111 2 1 
ai a 2 a 3 2a 4 5a5 5a6 
If we substitute the values (|3.77|l . (|3.79l) for a 3 , a 4 , a 2 , a6 into l|3.75|l we see that 
the coefficient of G(l) is bounded below by 

(3.82) 2.64- 1.6 G(2)- 5o (a 1 ), a = G(2)/G(l), 

where the function g a (z) is defined by 



, , 2 8a 
?a(z) = I + 5l5-4iI^ 



< z < 5/4, a > 0. 
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We can easily compute the minimum of g a to be 



inf g a (z) = 

0<z<5/4 



V5- 



/25. 



Using the fact that a < 1/5, G(2) < G(l)/5 < 1/10 we see that the quantity 
is bounded below by 2.64 - .16 - 2.304 > 0. 



□ 



Proof of Theorem 1.2. (d — 2, L\ = 4) : We need only verify that the function G 
defined by l|3.32[l . 1|3.53[) satisfies the inequalities (|3.57[1 . Since G(y), y > 1, decays 
exponentially one can verify these inequalities with aid of a computer. In particular 
we see that 

G(0) = .7071, G(l) = .1213, G(2) = .0208, 
correct to 4 decimal places, whence (|3.57|) holds. □ 

4. Proof of Theorem 1.2 

In this section we obtain the formula (|1.7|) of §1 for the effective diffusion constant 
which generalizes the formulas obtained in §3. We take L\ ~ 2L with L > 2 in 
Lemma 2.6. Then fl — {(n,y) : 1 < n < L, y e f2d-i}. For y £ Sld-i> 1 < j < L 
we define 5 j {y),5 j {y) by 



(4.1) 



We see from l|2.62|l . (I4.1|) that ip* satisfies the system of equations, 
(4.2) 

t*£V(l,y) - 6 2 (y)<p*(2,y) - 5i(i/)^(l >tf ) 



5i(y)v*(l,y) 



- 6 L (y)^*(L,y) 5 L _ 2 {y)y*{L - 2, y) 

- S L (y)v*(L,y) - 6 L -i(v)<P*(L-l,y) 



0, 
0, 
0, 
0, 
0, 



where A = A^-i is the d — 1 dimensional Laplacian. If we add all the equations in 
(|4.2I) we obtain the equation 



-Aj2ip*(j,y)=0, yeOi-i. 
On using the normalization (<P*)q = 1 we conclude that 

L 

(4.3) $>*Cj,») = A ye^-i- 

i=i 

Evidently we can rewrite the first equation of l|4.2|l as 
A 



(4.4) 



2d 



h(y) 



( p*(l,y)-5 2 (y) l p*(2,y) = 0. 



If we add (|4.4|l to the second equation of l|4.2|l we obtain the equation 



(4.5) 



A. 
2d 



2d 



•*2(y) 



^(2,y)-<5 3 (7/K(3,y) = 0. 
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Adding Q4.5p to the third equation of l|4.2p and proceeding similarly with subsequent 

equations we obtain the system, 

(4.6) 



A r / N 



A 



A 

Yd 



>L-1 



(y) 



<p*(L-l,y)-6 L (y)p*(L,y)=0, 



2dJ ^,y)+-+[- Yd )<P*(L-2,y) + 

where we have omitted the final equation of l|4.2|l . From 14. 411 . (I4.5|l . (|4.6I) we can 

write tp*(j,y), 2 < j < L, in terms of (p*(l,y). Substituting these into (|4.3|) we 
obtain an equation for <^*(1, y) of the form 

(4.7) £<p*(l,y)=L, t/GOd-i, 

where £ is an operator on functions on O^-i- 

Next we consider the equations (|2.48|) . i|2.53p for the function tfj on Cl. Thus ipo 
satisfies the system of equations, 
(4.1 

V>o(l,y) - 5i(y)^ (2,y) + 

^(y)V'o(3,y) - 5 2 (y)ip (l,y) 



(-A+2) , 



2d 



(^)^o(i-i,y) h-i(v)ML,v) 

t^+31 MLi y) + 6 L (y)ipa(L,y) 
We can rewrite the first equation of (|4.8|l as 

(-A + 4) 



S L -i(y)ML-2,y) 
Sl{v)iPo{L- l,y) 



0. 
0. 
0. 

o, 



(4.9) 



2d 



?/>o(l, y) = [Ml, y) + ^o(2, y)] = h{y)u{l, y), y G n d . 



where w(l) = ^o(l) + ^o(2). If we put now u(2) = "0o(3) — ^o(l) then on using 
(|4.9|l we see that the second equation of l|4. 8J) is the same as 

A -S 2 u{2) = 0. 



(4.10) 



2rf 



Observe that (|4. is identical to l|4.4p under the reflection b — > —6. We can 
similarly obtain the reflection of the equations (|4.5H , (|4.6|) by defining the variables 
u(j),j = 3, ... by 

= M.7 + i) - Mj ~ 1), J = 3, L—l, 

Let us assume that the u(j),j — 1, , J— 1, satisfy the reflection of the first J — 2 

of the equations (|4.4jl . Ij4.5|l . (|4.6|l . We show that u(J) then satisfies the (J — l)st 
equation provided J < L — 1. To see this we consider the Jth equation of (|4.8p 
which we may write as 

(4.11) { ^d 2 ) MJ) ~ ~ 5j [M(J) + MJ ~ 1)] ~ SjMJ ~ 1} = °' 

We may rewrite 114. lip as 



(4.12) 



A . 
-2d +6 ^ 



u{J - 1) - 6ju{J) + 8j-xu{J - 1)+ 
(-A + 2) 



2d 



2)--^ (J-l) = 0. 



HOMOGENIZATION OF A DIFFUSION PERTURBED BY A VECTOR FIELD 



33 



If J = 3 then wc have that 
5j-MJ-l) + 

5 2 u{2) + 

8 2 u{2) + 



(-A + 2) 2 

>-MJ-2)-— ^ (J-1) = 



2d 
(-A + 4) 
2d 



MV~Yd [Ml) + M2)] = 



8 x u{l)-—u{l)=5 2 u{2)-6 x u{l) 



2d 



■Si 



u(l) - S lU (l) 



2d 



u(l). 



We have shown that the result holds for J — 3. More generally we have that 



Sj-m{J-i) + 



(-A + 2) 



2d 



^o(J-2)-^ (J-l) = 



A 

2d 



(5,7-2 



u(J- 2) 



J-3 



2d 



+ E-^) + ^^^(J-2)-A^!./ -1) 
i=i 

• 7 - 2 a 2 
= E-^')-^-2 M ( J - 2 )-^o(J-3) 

(-A + 2) 



2d 



2d 



To complete the proof we need then to show that 

A,,^-^- A + 2 ) 

2d 



-6j- 2 u(J - 2) - --^ (J - 3) + v . ■' "V o(J - 2) = 0. 



2d 

This last equation is however simply the (J — 2)nd equation of 14.8(1 . 

We have shown that u(j),j = 1, L — 1 satisfies the reflection of the first L — 2 
equations of JOJ, 1)4. 5fl . (|4.6p . Define now u(X) = 1 - ^o(-k) - i>o(L - 1), whence 
there is the identity, 

L 

(4.13) Yl <l) = L 

We shall show that the u(j),j = 1, L satisfy the reflection of the final equation 
of (|4.t>p . To see this we write the final equation of 14.8fl as 

'-A + 2 N 



2d 



ML) +S L [1- u(L) - ML - 1)] - S L ML -1) = 6i 



whence we have that 
'-A + 2' 



2d 



[u(L - 1) + ML - 2)] - 6 L u(L) - —ML - 1) = 0. 



We may rewrite the previous equation as 



A . 
"2d + ^ 



u(L-l)-5 L u{L)+S L -iu(L-l)+ ( + 2) ML-2)-^ML~l) = 0. 



Now if we use the identity already established 
5 L -iu(L-l) = 



A . 
'2d +SL - 2 



(L-2) + L f-*uU), 
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we see that it is sufficient to show that 

8 L ^ 2 u(L - 2) + 2 ^ % o(L - 2) - ^{L - 1) = 0. 

This last equation is just the (L — 2)nd equation of (|4.8|) . 

Let Cr be the reflection of the operator C of l|4.7|) obtained by replacing b by 
—b. Then on comparing (|4.3|) . I|4.13[) we see that it(l) satisfies the equation 

(4.14) C R u(1) = 1. 

We are able now to come up with a new formula for the effective diffusion constant. 
On using (|2.54|) , l|4.7|) , l|4.9|) , l|4.14|l we have that the effective diffusion constant is 
given by 

(4.15) 8L 2 d < [SxC- 1 1] (-A + 4)- 1 [S^ 1 l] ) , 

where ( • ) is the uniform probability measure on fld-i- The formula l|1.7|) follows 
from (|4.15() . In order for (|4.15() to be valid we need to show that C is invertible. 

Lemma 4.1. Let C be the matrix defined by |^.7| ). Then C is invertible and the 
matrix C~ l has all positive entries. 

Proof. We proceed by induction. For k = 2, 3 let Ck be the operator (|4.7II when 

L = k. It is easy to see from (|4.r>[) - (14.611 that the Ck satisfy the recurrence relation, 



(4.16) C 



k+l 



A 

2d 



- <h 



k+l 



-C k -i, k>l; C o = 0, Ci 



J k+i 



The result will follow by showing that the matrices Ak = Ck-iCZ , k > 2, have all 
positive entries and principal eigenvalue strictly less than 1. Evidently this is the 
case for k — 2. Now from (|4.16|) we see that the Ak satisfy the recurrence relation, 



(4.17) 



Ai. + \ = { -— + S k + S k +i - 5kA k 



If Ak has all positive entries with principal eigenvalue strictly less than 1 then the 
matrix [—A/2d + 5k + 4+i] _ "kAk has the same property and the matrix Ak+i 
defined by (|4.17|) has all positive entries. To conclude the induction step we need 
therefore to show that Ak+i has principal eigenvalue strictly less than 1. To see 
this note that if 1 denotes the vector with all entries 1 then 



— + 8k + Sk+i - 8 k A k 



1 > 8 k+i, 



whence we conclude that 

A 

2d 



+ 8 k + S, 



k+l 



~8 k A k \ 4+i(l)<l. 



□ 
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